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Abstract 

Interwiners are the building blocks of spin-network states. The space of intertwiners is 
the quantization of a classical symplectic manifold introduced by Kapovich and Millson. 
Here we show that a theorem by Minkowski allows us to interpret generic configurations 
in this space as bounded convex polyhedra in R 3 : a polyhedron is uniquely described by 
the areas and normals to its faces. We provide a reconstruction of the geometry of the 
polyhedron: we give formulas for the edge lengths, the volume and the adjacency of its 
faces. At the quantum level, this correspondence allows us to identify an intcrtwincr with 
the state of a quantum polyhedron, thus generalizing the notion of quantum tetrahedron 
familiar in the loop quantum gravity literature. Moreover, coherent intertwiners result to 
be peaked on the classical geometry of polyhedra. We discuss the relevance of this result for 
loop quantum gravity. In particular, coherent spin-network states with nodes of arbitrary 
valence represent a collection of semiclassical polyhedra. Furthermore, we introduce an 
operator that measures the volume of a quantum polyhedron and examine its relation 
with the standard volume operator of loop quantum gravity. We also comment on the 
semiclassical limit of spinfoams with non-simplicial graphs. 
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1 Introduction 

Loop quantum gravity (LQG) is a continuous theory, whose Hilbert space is the direct sum 
of spaces associated to graphs T embedded in a three-dimensional hypersurface, % = ©r%r- 
It is often convenient to consider a single graph T, and the associated Hilbert space T~Lr- The 
truncation captures only a finite number of degrees of freedom of the theory. An important 
question for us is whether these degrees of freedom can be "packaged" as to provide some 
approximate description of smooth 3d geometries [1, 2]. We specifically think that it would be 
useful to have a picture of the classical degrees of freedom captured by Hr in terms of discrete 
geometries. Such knowledge is for instance relevant for the interpretation of semiclassical states 
restricted on "Hr- 

As it turns out, useful insights can be gained looking at the structure of It decomposes 
in terms of SU(2)-invariant spaces T~L f associated to each node of valence F. For a 4-valent node, 
it has been known for quite some time that an intertwiner represents the state of a "quantum 
tetrahedron" [3, 4], namely the quantization of the space of shapes of a flat tetrahedron in 
R 3 with fixed areas. For a generic valence F, a natural expectation would be a relation to 
polyhedra with F faces, as mentioned in [5] and [1]. In this paper we clarify the details of this 
correspondence. 

There are two keys to our result. The first one is the fact that Hp is the quantization of a 
certain classical phase space Sf, introduced by Kapovich and Millson in [6]. The second is the 
fact that there is a unique bounded convex polyhedron with F faces of given areas associated to 
each point of Sf- This is guaranteed by an old theorem by Minkowski [7]. The correspondence 
is up to a measure-zero subset of "degenerate" configurations, present also in the 4-valent case. 
Accordingly, we have the following relations: 

polyhedra with F faces < — > classical phase space Sf < — > intertwiner space Tip. 
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An immediate consequence of these results is a complete characterization of coherent states 
at a fixed graph: they uniquely define a collection of polyhedra associated to each node of 
the graph. This provides a simple and compelling picture of the degrees of freedom of 7ir 
in terms of discrete geometries, which are associated with a parametrization of the classical 
holonomy-flux variables in terms of the twisted geometries introduced in [1]. 

The paper is divided into two parts, concerning respectively the classical geometry of poly- 
hedra, and the notion of quantum polyhedron together with its relevance to loop gravity. The 
motivation for the first part comes from the fact that polyhedra have a rich classical geometry. 
One of the reasons why the notion of quantum tetrahedron has been so fruitful in the devel- 
opement of loop gravity and spinfoams is the fact that everybody understands the geometry 
of a classical tetrahedron. To make the extension to higher valence as fruitful, we need first of 
all to clarify a number of aspects of the geometry of polyhedra. 

Minkowski's theorem guarantees that a polyhedron can be reconstructed out of the areas 
and normals to its faces, just as it happens for the tetrahedron. The new feature here is 
that there are many possible polyhedra with the same number of faces which differ in their 
combinatorial structure, i.e. in the adjacency relations of the faces. In the first part of the 
paper (sections 2 and 3) we focus entirely on the classical geometry of polyhedra, and collect 
and in some cases adapt various results known in the mathematical literature. We discuss the 
combinatorial classes of polyhedra, and how the phase space of shapes at given areas can be 
divided into regions of different classes. We show explicitly how a given configuration of areas 
and normals can be used to reconstruct the polyhedron geometry, including its edge lengths, 
volume and combinatorial class. Furthermore, we discuss certain shape matching conditions 
which effectively restrict a collection of polyhedra to (a generalization of) Regge geometries. 

In the second part of the paper (sections 4 to 6) we discuss the quantum theory. We first 
review the construction of the quantization map between the phase space Sf and Hilbert space 
of intertwiners T-Lf- This leads to the interpretation of an intertwiner state as the state of 
a quantum polyhedron, and of coherent intertwiners [8, 9] as states describing semiclassical 
polyhedra. The relevance of polyhedra extend to the whole graph Hilbert space Hr, via the 
twisted geometries variables. The result provides an interpretation of coherent spin- network 
states in Hr as a collection of semiclassical polyhedra. 

Furthermore, we introduce a new operator which measures the volume of a quantum poly- 
hedron. Its definition is based on the knowledge of the classical system behind the intertwiner 
space T~Lf, and has the right semiclassical limit on nodes of any valence. We discuss its rela- 
tion with the standard volume operator of loop quantum gravity. Finally, we make some brief 
remarks on the polyhedral picture, Regge calculus and covariant spin foam models. 

2 The phase space of polyhedra 

2.1 Convex polyhedra and Minkowski theorem 

A convex polyhedron is the convex hull of a finite set of points in 3d Euclidean space. It can 
be represented as the intersection of finitely many half-spaces as 



V = { 



ni • x < hi, i = 1, . . . , m} 
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where rii are arbitrary vectors, and hi are real numbers. The abstract description (1) is non- 
unique and redundant: the minimal set of half-spaces needed to describe a polyhedron corre- 
sponds to taking their number m equal to the number of faces F of the polyhedron. In this 
paper we are interested in the description of a convex polyhedron with F faces in terms of vari- 
ables that have an immediate geometric interpretation: the areas of the faces of the polyhedron 
and the unit normals to the planes that support such faces. 

Let us consider a set of unit vectors rii G R 3 and a set of positive real numbers Ai such 
that they satisfy the closure condition 

F 

C = Y,Ani = 0. (2) 

i=l 

In the following, we will refer to this set as "closed normals". A convex polyhedron with F 
faces having areas Ai and normals rii can be obtained in the following way. For each vector 
rii consider the plane orthogonal it. Then translate this plane to a distance hi from the origin 
of R 3 . The intersection of the half-spaces bounded by the planes defines the polyhedron, 
rii ■ x < hi. We can then adjust the heights so that hi = hi(A) so that the faces have areas Ai. 

Remarkably, a convex polyhedron with such areas and normals always exists. Moreover, it 
is unique, up to rotations and translations. This result is established by the following theorem 
due to H. Minkowski [7, 10]: 

Theorem (Minkowski, 1897) 

(a) Ifni,...,riF are non-coplanar unit vectors and A±, . . . , Af are positive numbers such that 
the closure condition (2) holds, than there exists a convex polyhedron whose faces have 
outwards normals rii and areas Ai. 

(b) If each face of a convex polyhedron is equal in area to the corresponding face with parallel 
external normal of a second convex polyhedron and conversely, then the two polyhedra are 
congruent by translation. 

This unicity will play an important role in the following. Throughout the rest of the paper, we 
use simply polyhedra to refer to bounded convex polyhedra. 

2.2 Kapovich-Millson phase space as the space of shapes of polyhedra 

Let us consider F vectors in R 3 that have given norms A\ , . . . , Af and such that they sum up 
to zero. The space of such vectors modulo rotations has the structure of a symplectic manifold 
[6] and is known as the Kapovich-Millson phase space 1 Sf, 

S f = {ni G (S 2 ) F | Ef=i = 0}/SO(3) . (3) 

The Poisson structure on this 2(F— 3)-dimensional space is the one that descends via symplectic 
reduction from the natural SO(3)-invariant Poisson structure on each of the F spheres S 2 . 

1 In [6] it is also called the space of shapes of (bended) polygons. To be precise, it is a symplectic manifold 
up to a finite number of points, corresponding to configurations with one or more consecutive vectors collinear. 
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Action-angle variables for (3) are (F — 3) pairs with canonical Poisson brackets, 

{jUj, 9j} = 5ij. Here /ij is the length of the vector = Airai + . . . + Ai + irii + i (see Fig.l), and 
its conjugate variable Q\ is the angle between the plane identified by the vectors jli-i, fli and 
the plane identified by the vectors fli, fli+i- At fixed areas, the range of each fn is finite. 




Figure 1: A polygon with side vectors Ajrij and the (F— 3) independent diagonals. The space of possible 
polygons in R 3 up to rotations is a (2F — 6)-dimensional phase space, with action-angle variables the 
pairs of the diagonal lengths and dihedral angles. For non-coplanar normals, the same data 

defines also a unique polyhedron thanks to Minkowski's theorem. 

Thanks to Minkowski's theorem, a point in Sf with non-coplanar normals identifies a 
unique polyhedron. Accordingly, we refer to (3) as the space of shapes of polyhedra at fixed 
areas. Notice that (3) contains also configurations with coplanar normals: they can be thought 
of as "degenerate" polyhedra, obtained as limiting cases. The fact that the polyhedra with 
faces of given areas form a phase space will be important in section (4.1) where we discuss the 
Hilbert space of the quantum polyhedron. 

2.3 Classes of polyhedra with F faces 

The phase space Sf has a rich structure: as we vary the normals of a polyhedron keeping 
its areas fixed, not only the geometry, but in general also the combinatorial structure of the 
polyhedron changes; that is, the number of edges and the adjacency of faces. We refer to 
the combinatorial structure as the class of the polyhedron. In other words, there are two 
components to the shape of a polyhedron: its class, and its geometry (up to rotations) once the 
combinatorial structure is fixed. The different classes of polyhedra correspond to the different 
tessellations of a sphere having F faces. Which class is realized, depends on the specific value 
of the normals. This is a point we would like to stress: one is not free to choose a class, and 
then assign the data. It is on the contrary the choice of data that selects the class. This is 
an immediate consequence of Minkowski's theorem. Accordingly, the phase space Sf can be 
divided into regions corresponding to the different classes of polytopes with F faces. 

To visualize the class of a polyhedron it is convenient to use Schlegel diagrams [11, 10]. The 
Schlegel diagram of a polyhedron is a planar graph obtained choosing a face /, and projecting 
all the other faces on / as viewed from above. See Fig. 2 for examples. 

To understand the division of Sf into regions of different class, let us first give some 
examples, and postpone general comments to the end of the Section. In the most familiar 
F = 4 case, there is no partitioning of £4: there is a unique tessellation of the sphere, the 
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Figure 2: Some examples of Schlegel diagrams. From left to right, a tetrahedron, a pyramid, a cube 
and a dodecahedron. 

tetrahedron, and it is well known that there is always a unique tetrahedron associated with 
four closed normals. The first non-trivial case is F = 5, where there are two possible classes: a 
triangular prism, and a pyramid (see Fig. 3). Consider then the phase space S§. Minkowski's 



Dominant: 




Figure 3: Polyhedra with 5 faces: the two possible classes are the triangular prism (left panel) and the 
pyramid (right panel). The two classes differ in the polygonal faces and in the number of vertices. 

theorem guarantees that the same set (Ai,rii) cannot be associated to both classes, thus each 
point in S§ corresponds to a unique class. One might at first think that S§ can be more or less 
equally divided among the two classes, but this is not the case. In fact, notice that the pyramid 
is just a special case of the prism, obtained by collapsing to a point one of the edges connecting 
two triangular faces. The existence of a pyramid then requires a non-trivial condition, i.e. the 
presence of a 4-valent vertex. A moment of reflection shows that this condition can be imposed 
via an algebraic equation on the variables. Hence the shapes corresponding to pyramids span a 
codimension one surface in £5. Generic configurations of areas and normals describe triangular 
prisms, and the pyramids are measure zero special cases. We call dominant the class of maximal 
dimensionality, e.g. the triangular prism here. 

Let us move to F = 6, a case of particular interest since regular graphs in R 3 are six-valent. 
There are seven different classes of polyhedra, see Fig. 4. The most familiar one is the cuboid 
(top left of Fig. 4), with its six quadrilateral faces. Remarkably, there is a further dominant 
class: it is a "pentagonal wedge", i.e. a polyhedron with two triangles, two quadrilaterals and 
two pentagons as faces (to visualize it, immagine a triangular prism planed down on a corner, 
so that a vertex is replaced by a triangle). The remaining five classes are subdominant, because 
non-trivial conditions are required for their existence. Subdominant classes have fewer vertices 
and thus can be seen as special cases with certain edges of zero length. 2 

2 Among these, notice the class of codimension 3. It has six triangular faces and three four-valent vertices. 
This class is interesting in that it can be seen as two tetrahedra glued along a common triangle. Two arbitrary 
tetrahedra are denned by 12 independent numbers. In order for them to glue consistently and generate this 
polyhedron, the shape of the shared triangle has to match. This shape matching requires three conditions (for 
instance matching of the edge lengths), thus we obtain a 9-dimensional space of shapes. For fixed external areas, 
this is precisely the codimension 3 subspace in 56- Hence this class is a special case of two tetrahedra where 



6 



Dominant: 
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Figure 4: The seven classes of polyhcdra with 6 faces, grouped according to the dimensionality of their 
configurations. 

From the above analysis, we expect that the phase space Sq can then be divided into 
regions corresponding to the two dominant classes, separed by the subdominant ones. This 
is qualitatively illustrated in Figure 5. To confirm this picture, we performed some numerical 




Figure 5: Pictorial representation of the phase space: it can be mapped into regions corresponding to 
the various dominant classes (two in the example). The subdominant classes separe the dominant ones 
and span measure-zero subspaces. 

investigations. Using the reconstruction algorithm, which we introduce in the next Section, we 
can assign a class to each point in Sq. In Fig. 6 we give an explicit example of a 2d and a 3d 
slice of the 6d space Sq, which shows the subdivision into the two dominant classes. 

After this brief survey of some specific examples, let us make some general statements. 

• The phase space Sf can be divided into regions corresponding to different classes. The 
dominant classes, generically more than one, cover it densely, whereas the subdominant 
ones span measure-zero subspaces. The dominant classes in phase space correspond to 
polyhedra with all vertices three- valent, that is the dual to the tessellation is a triangu- 
lation. This condition maximizes both the number of vertices, V = 3(F — 2), and edges, 
E = 2{F — 2). Subdominant classes are special configurations with some edges of zero 
lengths and thus fewer vertices. 

• Since all classes correspond to tessellations of the sphere with F faces, they are connected 
by Pachner moves [12]. The reader can easily find a sequence of moves connecting all 

conditions are imposed for them to glue consistently. 
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Figure 6: Mappings of subspaces of 56 realized using the reconstruction algorithm of Section 3 and 
Wolfram's Mathcmatica. We subdivided the phase space into a regular grid, and had Mathematica 
computing the adjacency matrix of the area-normal configurations lying at the center of the cells. This 
associates a unique class to each cell of the phase space. The information is colour-coded, cuboids in 
blue, pentagonal wedges in red. With this mapping of finite resolution we have measure-zero probability 
of hitting a subdominant class, thus the latter are absent in the figures. The holes arc configurations 
for which our numerical algorithm failed. Concerning the specific values of the example, the areas are 
taken to be (9, 10, 11, 12, 13, 13). In the left panel, we fixed /ii = 15, Q\ = y^7r, ^t 2 = 13, 2 = j^it, and 
plotted the remaining pair (1^3,63). In the right panel, we fixed \ii = (15, 13, 17) and plotted the three 
angles 0j. 



seven classes of Fig. 4. To start, apply a 2-2 move to the upper edge of the inner square 
of the cuboid to obtain the pentagonal wedge. 

• The lowest-dimensional class corresponds to a maximal number of triangular faces, a 
condition which minimizes the number of vertices. When all the faces are triangular, the 
polyhedron can be seen as a collection of tetrahedra glued together, and with matching 
conditions imposed along all shared internal triangles. 

2.3.1 Large F and the hexagonal dominance 

The number of classes grows very fast with F (see for instance [13] for a tabulation). In the 
examples above with small F, we have been able to characterize the class looking just at how 
many faces have a certain valence. However as we increase F we find classes with the same 
valence distribution, but which differ in the way the faces are connected. To distinguish the 
classes one needs to identify the complete combinatorial structure of the polyhedron. This 
information is captured by the adjacency matrix, which codes the connectivity of the faces of 
the polyhedron. Below in Section 3.3 we will show how this matrix can be explicitly built as a 
function of areas and normals, and give some explicit examples. 

An interesting question concerns the average valence of a face, defined as {p) = 2E/F. A 
simple estimate can be given using the fact that the boundary of any polyhedron is a tessellation 
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of the two-sphere, therefore by the Euler formula F — E + V = 2. For the dominant classes, 
which are dual to triangulations, the additional relation 2E = 3V holds, hence E = 3(F — 2) 
and we get [p) = 6(1 — 2/F). For large F, we expect the polyhedron to be dominated by 
hexagonal faces. This expectation is immediately confirmed by a simple numerical experiment. 
The specimen in Figure 7, for instance, has F = 100 and (p) ~ 5.88. Notice also from the 



Figure 7: A polyhedron with F = 100 drawn with Wolfram's Mathcmatica, using the reconstruction 
algorithm. The example has all areas equals and normals uniformly distributed on a sphere. Notice 
that most faces have valence 6, and that triangles are nowhere to be seen. 

image that there are no triangular faces, consistently with the fact that they tend to minimize 
the number of vertices and are thus highly non-generic configurations. 

3 Polyhedra from areas and normals: reconstruction procedure 

So far we have discussed how a point in Sf specifies a unique polyhedron, and the existence 
of different combinatorial structures. We now describe how the polyhedron can be explicitly 
reconstructed from areas and normals. The reconstruction will allow us to evaluate completely 
its geometry, including the lengths of the edges and the volume, and to identify its class through 
the adjacency matrix, thus being able to associate a class with each point of Sf- 

The main difficulty in developing a reconstruction algorithm is that, given the areas and 
the normals, it is not known a priori which faces of the polyhedron are adjacent. The adjacency 
relations of the faces (and the combinatorial class of the polyhedron) are to be derived together 
with its geometry. This can be done in two steps. The first step uses an algorithm due to 
Lasserre [14] that permits to algebraically compute the lengths lij(h,n) of all the edges of the 
polyhedron as defined by hi and n«, as in (1). The second step consists of solving a certain 
quadratic system to obtain the values of the heights hi for given areas. 

3.1 Lasserre's reconstruction algorithm 

We now review Lasserre's procedure, and adapt it to the three-dimensional case of interest 
here. The basic idea of the reconstruction algorithm is to compute the length of an edge as 
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the length of an interval in coordinates adapted to the edge. Consider the i-th face. From the 
defining inequalities (1), we know that points i£l 3 on this face satisfy 

rii ■ x = hi (4a) 
rij • x < hj, i ^ j. (4b) 

We consider the generic case in which rii • rij ^ ±1 Vi, j (these special configurations can be 
obtained as limiting cases). We introduce coordinates fji adapted to the face, that is 

m-yi = 0, yi = x - (x ■ m)ni. (5) 

Using (4a) we get x = hirii + m, which inserted in (4b) gives 

yi-nj<rij, i^j, (6) 

where we have defined 

rij = hj - (m ■ rij)hi . (7) 
Hence, the i-th face can be characterized either in terms of the x or the yi coordinates, 

x ■ rii = hi \ yi ■ Hi = 



(8) 

rij ■ x < hj, i^j \ yi-rij < rij{h, n), i / j 



Notice that r^j \J\ — (rii ■ rij) 2 is the distance of the edge ij from the projection of the origin 
on the i-th. face. 

The next step is to iterate this process and describe an edge in terms of its adapted coordi- 
nates. We start from the i-th. face again, and assume that it is connected to the face j, so that 
the two faces share an edge. Points on the edge ij between the i-th and the j-th face satisfy 

yi ■ rii = (9) 

Vi - n 3 = r ij ( 10 ) 

yi-n k < r ik , k / i,j. (11) 
As before, we introduce coordinates Zij, adapted to the edge, 

m ■ = rij ■ = 0, z^ =yi- [rij - (m ■ nj)rii] —f Uj . (12) 

1 (rii ' Hj) 

Using (10) we get that for a point in the edge 

y% = [nj - {rii ■ nj)m\ -, ^r + z ij- ( 13 J 

1 — (rii • nj) 

Plugging this in (11) gives 

Zij ■ n k < b ij>k , (14) 

where we have defined 

7 u t \u {nj ■ n k ) - (rii ■ nj)(m ■ n k ) u( v, MtC . 

hj,k = h k - (rii ■ n k )hi J - [hj - h^rii ■ n,-) . (15) 

1 — (rii • n j) 
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Summarizing as before, going to adapted coordinates the edge is defined by 
y i -n i = ( -ni = 

Vi-rij =nj(h,n) — > < Zij-rij = (16) 
Vi ■ n k < r ik (h,n), k^i,j. [ z tj ■ n k < b ij)k (h,n), i^j^k 

At this point we are ready to evaluate the length of each edge. To that end, we parametrize 
the coordinate vector in terms of its norm, say A, and its direction which is given by the 
wedge product of the two normals, 

Zij = A , (17) 
y/1 - (rii ■ rij) 2 

If we define 

niAnj-rik 

<Hj,k = 7j 7 W > ( 18 ) 

a/1 - (rii • rij) 2 

we can rewrite the inequalities in (16) as 

Aa ijjfc < b ij)k . (19) 

Finally, the length of the edge is the length of the interval determined by the tightest set 
of inequalities, i.e. 

min [hhh\- max ihhhX. (20) 



fc|iij,fc>0 {Clij yk ) k\a ijtk <0 {Clij jk 

Here the minimum is taken over all the fc's such that clm k is positive, and the maximum over all 
the A;'s such that is negative. This quantity is symmetric [14] and satisfies a key property: 
it can be defined for any pair of faces ij, not only if their intersection defines an edge in the 
boundary of the polyhedron, and it is negative every time the edge does not belong to the 
polyhedron [14]. Thanks to this property, we can consistently define the edge lengths for any 
pair of faces ij as 

t%j(h, n) = max i 0, min < — 1 — max | \ 1 . (21) 

The result is a matrix whose entries are the edge lengths (as a functions of the normals and the 
heights) if the intersection is part of the boundary of the polyhedron, and zero if the intersection 
is outside the polyhedron. 

This formula completes Lasserre's algorithm, and permits one to reconstruct the polyhedron 
from the set (hi, rii). To achieve a description in terms of areas and normals, we need one more 
step, that is an expression for the heights in terms of the areas. This can be done using (21) to 
compute the areas of the faces. We consider the projection of the origin on the face, and use 
it to divide the face into triangles. Recall the Lasserre's procedure has provided us with the 
distance between an edge and the projected origin, see (8). We thus can write 

1 A r- 

Ai = -Y , 13 in . (22) 
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Notice that both rij(h,n) from (7) and £ij(h,n) from (21) are linear in the heights. Hence, the 
area is a quadratic function, 

F 

Ai(h, n)=J2 M t( n i> ■ • • > ri F )hjh k , (23) 

j,k=l 

where Mj is a matrix depending only on the normals. This homogeneous quadratic system can 
be solved for hi(A, n). The existence of a solution with hi > Vz is guaranteed by Minkowski's 
theorem. However, the solution is not unique: in fact, we have the freedom of moving the 
origin around inside the polyhedron, thus changing the value of the heights without changing 
the shape of the polyhedron. A method which we found convenient to use is to determine a 
solution minimizing the function 

/(h^^Y^iMK^-Aif (24) 

i 

at areas and normals fixed, with Ai{h,n) given by (23). This is the method used in the 
numerical investigations of Figs. 6 and 7. 3 

Finally, from the inverse we derive the lengths as functions of areas and normals, which 
with a slight abuse of notation we still denote in the same way, 

e ij (A,n)=£ ij (h(A,n),n). (25) 

These expressions are well-defined and can be computed explicitly. 



3.2 Volume of a polyhedron in terms of areas and normals 

Let us call V(Ai,rii) the convex subset of M 3 corresponding to the polyhedron. Its volume is 
simply the integral on this region of the Euclidean volume density: 

V{Ai,m)= [ d 3 x. (26) 
Jv(Ai,m) 

An interesting question is how to compute efficiently the volume integral (26). The simplest 
way is to use the algorithm described in the previous section: we chop the region V(Ai, rtj) into 
pyramids with a common vertex in its interior and bases given by the faces of the polyhedron. 
In this way the volume is just the sum of the volumes of the pyramids, i.e. 

i f 

V{A i ,n i ) = -Y J h i A i . (27) 
6 i=i 

Here hi = hi(A, n) are the heights of the pyramids expressed in terms of the areas and normals 
via Lasserre's algorithm. 

3 Concerning Fig. 6, we can also give now more details on the holes: these are configurations for which 
the numerical algorithm to solve (23) failed. This limitation can be easily improved with a better inversion 
algorithm, or by choosing a configuration slighly off the center of the cell. 
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The volume can be used to define a volume function on the phase space Sf- To that end, 
notice that (27) is not defined for configurations with coplanar normals, which on the other 
hand do enter Sf- However, it can be straightforwardly extended to a function on the whole 
Sf by defining it to be zero for coplanar configurations. Furthermore, the resulting phase space 
function is continuous. 4 Since the volume is manifestly invariant under rotations, it can also 
be written as a function of the reduced phase space variables only, that is, V(Ai, fif~, 0/.). To 
do so explicitly, one uses the relation m = ni which is straightforward to derive once a 

reference frame is chosen. 

The volume of the polyhedron as a function of areas and normals has a number of interesting 
properties: 

CI. Non-negative phase-space function. The volume is by construction non-negative, and at 
given areas, it vanishes only when the normals rij lie in a plane. This in particular implies 
that the volume vanishes for F = 2 and 3. 

C2. Boundedness. For fixed areas Ai, the volume is a bounded function of the normals. We 
call V max (Ai) the volume of the polyhedron with maximum volume, 5 

V max (A { ) = sup{ V (A { , m ) } . (28) 

In particular, V max (Ai) is smaller that the volume of the sphere that has the same surface 
area as the polyhedron. Therefore we have the bound 

o<V(Ai,m)< ( ^^ 2 • (29) 

3V47T 

C3. Face- consistency. If we set to zero one of the areas such that the result is still a non- 
degenerate polyhedron, the function (27) automatically measures the volume of the re- 
duced polyhedron with F — 1 faces. 

In conclusion, a point in Sf determines uniquely the whole geometry of a polyhedron and 
in particular its edge-lengths £ij (21) and its volume (27). 6 Now we show how these data can 
be used to identify the class of the polyhedron. 

4 In order to to see this, one shows that the limit of coplanar normals exists and the volume tends to zero in 
this limit. From property (C3) - see below, a general F-valent coplanar configuration can be obtained from a 
F + 1 configuration in the limit of zero base's area. 

5 Notice that there can be more than one polyhedron that attains maximum volume. For instance, in the case 
F = 4, there are two parity-related tetrahedra with maximal volume. 

6 It is worth adding that the problem of computing the volume of a given polyhedron is a complex and well 
studied topic in computational mathematics [15, 16], hence better procedures than the one used here could in 
principle be found. However, the usual starting point for common algorithms is the knowledge of the coordinates 
of vertices, or the system of inequalities (1). Therefore the methods need to be adapted to obtain formulas in 
terms of areas and normals. The main difficulty is clearly that the adjacency relations of the faces are to 
be derived together with the geometry. We found Lasserre's algorithm to be the most compatible with these 
necessities, thanks to the fact that the lengths are reconstructed algebraically. Numerical algorithms for the 
volume and shape reconstruction from areas and normals are developed in the study of extended Gaussian 
images in informatics [17], however there are no analytical results. 
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3.3 Adjacency matrix and the class of the polyhedron 

The adjacency matrix A of the polyhedron is defined as 



Ai 



1 if the faces i and j are adjacent 
otherwise 



hJ 



1,...,F 



(30) 



Notice that Aij coincides with the matrix l^j in (21) with all the non-zero entries normalized 
to 1: the recontruction algorithm gives us the adjacency matrix for free. 

The symmetric matrix Aij contains information on the connectivity of the faces as well as 
on the valence of each face, thus the class of the polyhedron can be identified uniquely from it. 
The valence pi of the face i can be extracted taking the sum of the columns for each row, 



Pi — ^2 Aij . 

3=1 

For example, for the two classes with F = 5 of Fig. 3 we have 



(31) 




A 



( 1 1 1 \ 
111 
110 11 
1110 1 

V i i i i o ) 



p= (3,3,4,4,4) 




A 



( 



\ 



\ 



J 



p= (4,3,3,3,3) 



From graph theory [18], we known that (30) has a number of interesting properties that 
can be related to the geometrical parameters of the polyhedron. For instance, the number of 
walks from the face i to the face j of length r is given by the matrix elements of the r-th power 
(A r )ij. From this property we deduce that the number E of edges of the polyhedron is 



E 



1 



TiA 



(32) 



This expression generalizes the value E = 3(F — 2) valid for the dominant classes. 

Higher traces are related to the number of loops of a given lengh. For instance, the number 
of closed loops of length 3 is given by (1/6) Tr^ 3 . 

Through the adjacency matrix, obtained via the reconstruction procedure, areas and nor- 
mals identify a unique class, and thus permits the division of Sf. 
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3.4 Shape-matching conditions 



Knowing the complete geometry of the polyhedra allows us also to address the following sit- 
uation. Suppose that we are given two polyhedra in terms of their areas and normals, and 
that we want to glue them by a common face. Even if we choose the area of the common face 
to be the same, there is no guarantee that the shape of the face will match: The two sets of 
data will in general induce different shapes of the face. That is, the face has the same area but 
it can be two different polygons altogether. In order to glue the polyhedra nicely, one needs 
shape matching conditions guaranteeing that the shared face has the same geometry in both 
polyhedra. 

If both polyhedra are tetrahedra, the problem has been solved in [19] . One uses the fact that 
the shape of the common triangle matches if two lengths, or two internal angles, are the same. 
The internal angles a can be expressed in terms of the 3d dihedral angles of the tetrahedron 
as follows, 



cos a 



jk 



COS + COS (pik COS 4>jk 

sin cf> ik sin 4> jk 



(33) 



Here the faces i, j and k all share a vertex, and a l - k is the angle between the edge ij and the 
edge ik inside the triangle i. Consider now the adjacent tetrahedron. Its geometry induces for 
the same angle the value 



COS (Xji k i 



cos •/ + cos i 



4' cos M 



sm< 



ik 



, sin i 



jk' 



Hence, for the shape to match it is sufficient to require 



j'k' 



(34) 



Ckl,ij{4>) 



cos a 



cos a 



j'k' 







(35) 



for two of the three angles of the triangle. These shape matching conditions are conditions 
on the normals of the two tetrahedra. See left panel of Figure 8 for an illustration of these 
relations. 





Figure 8: The geometric meaning of equation (35): the 2d angle ay,^ belonging to the shaded triangle 
can be expressed in terms of 3d angles associated the thick edges of the tetrahedron k, or equivalenty 
of the tetrahedron /. 
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The simplicity of the conditions (35) is a consequence of the fact that two triangles with 
the same area are congruent if two angles match. For the general case, the face to glue is now 
a polygon and the number of conditions greater. One needs to make sure that the valence p 
of the polygon is the same. Then, the number of independent parameters of a polygon on the 
plane is 2p — 3, hence giving the edge lengths is not enough, and p — 2 additional conditions 
are needed. A convenient procedure is the following. Identify the faces of the two polyhedra 
that, having the same area, we want to match. From the reconstruction algorithm, we know 
the edge lengths £y of the face viewed from one polyhedron. Then, for all j such that £ij ^ 0, 
we consider the face normals rij projected on the plane of the i-th face, 

rij - (rii ■ nArii rij - cos (fon^ 

nj = 7- 1 - f Vt = ■ a ■ ( 36 ) 

\7ij — (n» • nj)ni\ smcpij 

The set (£j,hj) defines a unique polygon in the plane identified by rij, thanks to a two- 
dimensional version of Minkowski's theorem. Then, we do the same with the second poly- 
hedron, obtaining a second set (£'j,n'j) living in the plane identified by n\. Finally, the shape 
matching conditions consist of imposing the equivalence of these two flat polygons up to rota- 
tions in three-dimensional space. Notice that the shape matching are now conditions on both 
the normals and the areas of the two polyhedra. 



4 Relation to loop quantum gravity 

Thus far we have been discussing classical properties of polyhedra. In the rest of the paper, 
we discuss the relevance of polyhedra for loop quantum gravity. The relation comes from the 
following two key results: 

(i) Intertwiners are the building blocks of spin-network states, an orthonormal basis of the 
Hilbert space of loop quantum gravity [20, 21] 

(it) Intertwiners are the quantization of the phase space of Kapovich and Millson [22, 9, 23] 
(see also [24, 25]), i.e. of the space of shapes of polyhedra with fixed areas discussed in 
the previous sections. 

Therefore an intertwiner can be understood as the state of a quantum polyhedron, and spin- 
network states as a collection of quantum polyhedra associated with each vertex. 

In this section we review how (ii) and the notion of quantum polyhedron are established, 
observe that coherent intertwiners are peaked on the geometry of a classical polyhedron and 
discuss the relevance of this fact for the relation between semiclassical states of loop quantum 
gravity and twisted geometries. 



4.1 The quantum polyhedron 

Let us consider the space of vectors in 3d Euclidean space with norm j. This is a phase space, 
the Poisson structure being the rotationally invariant one proper of the 2-sphere of radius j. 



16 



As is well known, its quantization 7 is the representation space of SU(2). We are interested 
in the phase space Sf, that is the space of F vectors that sum to zero, up to rotations. The 
Poisson structure on Sf is obtained via the symplectic reduction of the Poisson structure 
on the product of F spheres of given radius. Thanks to Guillemin-Sternberg's theorem that 
quantization commutes with reduction, 8 we can quantize first the unconstrained phase space 
XiSj., and then reducing it at the quantum level extracting the subspace of ®iV"*' that is 
invariant under rotations. This gives precisely the intertwiner space Hp = I nv [®fLi V^ 3 '"]. 
The situation is summarized by the commutativity of the following diagram, 

XiSl — > ®iV* 

Symplectic reduction J, \, Quantum reduction 

S F — »■ U F 

The correspondence between classical quantities and their quantization is the following: up 
to a dimensionful constant, the generators Jj of SU(2) acting on each representation space 
are understood as the quantization of the vectors Ain%. In LQG the dimensionful constant is 
chosen to be the Immirzi parameter 7 times Planck's area SirLp, 

Aim — ► Ei = 8vr 7 Lp X. (37) 
The closure condition (2) on the normals of the polyhedron is promoted to an operator equation, 

F 

^2^=0. (38) 
i=i 

This condition defines the space of intertwiners, and corresponds to the Gauss constraint of 
classical General Relativity in Ashtekar-Barbero variables. 

One can then proceed to associate operators to geometric observables through the quanti- 
zation map (37). The area of a face of the quantum polyhedron is 



A i = \jE i -E i = toryL 2 P y/3i{ji + l) (39) 

and produces an equispaced quantization of the area Ai ~ ji for large spins, i.e. up to quantum 
corrections. Notice that an ordering can be chosen so that the area is exactly Ai = 8Tvy L p ji. 
This ordering will be considered below to simplify the construction of the volume operator. 

The scalar product between the generators of SU (2) associated to two faces of the polyhe- 
dron measures the angle 9ij between them [27], 

On = arccos — = =j . (40) 

3 v'./H.y, • lUM./, • 1) 



7 Notice that, as usual, the quantum theory requirers the quantization of some classical quantities. In this 
case the norm of the vector has to be a half-integer j, the spin. 

8 For the general theory see [26], for details on the application to the current system see [4] and in particular [9]. 
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Notice that the angle operators do not commute among themselves, therefore it is not possible 
to find a state for a quantum polyhedron that has a definite value of all the angles between its 
faces. Moreover, the adjacency relations of the faces is not prescribed a priori, thus B%j might not 
even be a true dihedral angle of the polygon. Therefore an eigenstate of a maximal commuting 
set of angles is far from the state of a classical polyhedron: it is an infinite superposition of 
polyhedra of different shapes, including different combinatorial classes. Semiclassical states for 
a quantum polyhedron are discussed in the next section. 

4.2 Coherent intertwiners and semiclassical polyhedra 

Coherent intertwiners for Hp were introduced in [8] and furtherly developed in [9, 23] (for 
previous related work, see [28]). These Livine-Speziale (LS) coherent intertwiners are defined 
as the SU(2)-invariant projection of a tensor product of states \ji,n>i) € 



The states \j,n) are SU(2) coherent states peaked on the direction n of the spin [29, 30], 



In (41), the unit-vectors n.j can be assumed to close, Yliji n i = 0- The reduced states are still 
an overcomplete basis of Hf, as a consequence of the Guillemin-Sternberg theorem [9, 31]. 

Coherent intertwiners are semiclassical states for a quantum polyhedron: the areas are 
sharp, and the expectation value of the non-commuting angle operators 9ij reproduces the 
classical angles between faces of the polyhedron in the large spin limit, 



Moreover, the dispersions are small compared to the expectation values. 

A useful fact is that coherent intertwiners can be labeled directly by a point in the phase 
space Sf of Kapovich and Millson, and therefore by a unique polyhedron. This provides a 
resolution of the identity in intertwiner space as an integral on Sf- To realize this reduction, 
it is convenient to parametrize Sf via F — 3 complex numbers Zk, instead of Let us 

choose an orientation in R 3 and consider the stereographic projection z% of the unit-vectors rij 
into the complex plane. 9 The F — 3 complex variables Zk are the cross-ratios [9] 




(41) 




(42) 




(43) 



Z k 




k 




(44) 



9 



The relation between the unit- vector n = (n x , n v ,n z ) and the stereographic projection is 



n x — in 



'4 



z = — 



1 — n~ 



= — tan - e 



where 8 and <j> are the zenith and azimuth angles of S 2 , and we have chosen to project from the south pole. 
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Given an orientation in R 3 , a set of normals rtj that satisfy the closure condition (2) can be 
obtained function of the cross-ratios, 

m = ni(Z k ) ■ (45) 

Coherent intertwiners can then be obtained via geometric quantization [9]: they are labeled 
by the variables Z k , that is \ji,Z k ), and are equal to the states ||jj,nj)| = \ji,n>i(Zk)) up to 
a normalization and phase. 10 The resolution of the identity is given by an integral over the 
variables Z k , 

^H F = d/x(Z fe ) \ji,Z k )(ji,Z k \ , (46) 

where the integration measure d/j,(Zk) = Kj i (Z k , Z k )Y\ k & 2 Zk depends parametrically on the 
spins ji and is given explicitly in [9] . The relevance of this formula for the following discussion 
is that it provides a resolution of the identity in intertwiner space as a sum over semiclassical 
states, each one representing a classical polyhedron: the intertwiner space can be fully described 
in terms of polyhedra. 11 



4.3 Coherent states on a fixed graph and twisted geometries 

The states \ji, Z k ) provide coherent states for the space of intertwiners only, and should not be 
confused with coherent spin-network states for loop quantum gravity. Nevertheless, classical 
polyhedra and coherent intertwiners are relevant to the full theory, as we now discuss. 

To relate polyhedra to loop quantum gravity, consider a truncation of the theory to a 
single graph T, with L links and iV nodes. The associated gauge-invariant Hilbert space %y = 
L2[SU(2) L /SU(2) Ar ] decomposes in terms of intertwiner spaces 7iF(n) = Inv^ienVW] as 

n r = ® k (8 n ^ (n) ) • (48) 

This Hilbert space is the quantization of a classical space 12 Sr = T*SU(2) i //SU(2) Ar , which 
corresponds to (gauge-invariant) holonomies and fluxes associated with links and dual faces of 
the graph. The double quotient // means symplectic reduction. The key result is that this 
space admits a decomposition analogous to (48). In fact, it can be parametrized as the following 
Cartesian product [1], 

Sr = X T * S * X S F(n), (49) 

l n 

10 The states \ Zk) also define an holomorphic representation of the quantum algebra of functions ip{Zk) = 
(ji,Zk\ip), see [23]. We will not use this representation in this paper. 

Recently [5, 32, 33] attention has been given to a second space for which polyhedra are relevant. This is a 
sum of intertwiner spaces such that the total spin is fixed, 

Uj= Inv[®? =1 V Ui) ]. (47) 

31 ■ -3F 

T.iH = -' 

The interest in this space is that it is a representation of the unitary group U(F). Vectors in this space represent 
quantum polyhedra with fixed number of faces and fixed total area, but fuzzy individual areas as well as shapes 
as before. Coherent states for (47) can be built using U(F) coherent states [32]. These are also peaked on 
classical polyhedra like the LS states (41), thus the results in this paper are relevant for them as well. 
12 Again, this is a symplectic manifold up to singular points [34]. 
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where T*S is the cotangent bundle to a circle, F the valence of the node n, and Sf is the 
phase space of Kapovich and Millson. 

The parametrization is achieved through an isomorphism between holonomy-fluxes and a 
set of variables dubbed "twisted geometries" . These are the assignment of an area A\ and an 
angle £j to each link, and of F normals rii, satisfying the closure condition (2), to each node. See 
[1, 2] for details and discussions. In this parametrization, a point of Sr describes a collection 
of polyhedra associated to each node. The two polyhedra belonging to nodes connected by a 
link I share a face. The area of this face is uniquely assigned to both polyhedra A\ (notice that 
this fact alone does not imply that the shape of the face matches - more on this below). The 
extra angles £j carry information on the extrisic geometry between the polyhedra. 

The isomorphism (49) and the unique correspondence between closed normals and polyhe- 
dra means that each classical holonomy-flux configuration on a fixed graph can be visualized as 
a collection of polyhedra, together with a notion of parallel transport between them. Just as the 
intertwiners are the building blocks of the quantum geometry of spin networks, polyhedra are 
the building blocks of the classical phase space (49) in the twisted geometries parametrization. 

What is the relevance of this geometric construction to the quantum theory? Coherent 
states for loop quantum gravity have been introduced and extensively studied by Thiemann 
and collaborators [35, 36, 34]. Although the states for the full theory have components on each 
graph, one needs to cut off the number of graphs to make them normalizable. In practice, it 
is often convenient to truncate the theory to a single graph. This truncation provides a useful 
computational tool, to be compared to a perturbative expansion, and has found many appli- 
cations, from the study of propagators [37] to cosmology [38]. In many of these applications, 
control of the semiclassical limit requires a notion of semiclassical states in the truncated space 
7ir- The truncation can only capture a finite number of degrees of freedom, thus coherent 
states in Hr are not peaked on a smooth classical geometry. Twisted geometries offer a way 
to see them as peaked on a discrete geometry, to be viewed as an approximation of a smooth 
geometry on a cellular decomposition dual to the graph T. The above results provide a com- 
pelling picture of these twisted geometries in terms of polyhedra, and thus of coherent states 
as a collection of semiclassical polyhedra. 

There is one subtlety with this geometric picture that should be kept in mind, which justifies 
the name "twisted" geometries: they define a metric which is locally flat, but discontinuous. To 
understand this point, consider the link shared by two nodes. Its dual face has area proportional 
to A\. However, the shape of the face is determined independently by the data around each 
node (i.e. the normals and the other areas), thus generic configurations will give two different 
shapes. In other words, the reconstruction of two polyhedra from holonomies and fluxes does 
not guarantee that the shapes of shared faces match. Hence, the metric of twisted geometries 
is discontinuous across the face [1, 2]. 13 See left panel of Figure 8. 

One can also consider a special set of configurations for which the shapes match, see right 
panel of Figure 8. This is a subset of the phase space Sr where the shape matching conditions, 
discussed earlier in Section 3.4, hold. This subset corresponds to piecewise flat and continuous 
metrics. For the special case in which all the polyhedra are tetrahedra, this is the set-up 

13 Aspects of this discontinuity have been discussed also in [39, 40] 
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of Regge calculus, and those holonomies and fluxes indeed describe a 3d Regge geometry: 
twisted geometries with matching conditions amount to edge lengths and extrinsic curvature 
dihedral angles [1, 2]. This relation between twisted geometry and Regge calculus implies that 
holonomies and fluxes carry more information than the space of Regge calculus. This is not in 
contradiction with the fact that the Regge variables and the LQG variables on a fixed graph 
both provide a truncation of general relativity: simply, they define two distinct truncations of 
the full theory. See [2] for a discussion of these aspects. 

For an arbitrary graph, the shape-matching subset describes a generalization of 3d Regge 
geometry to arbitrary cellular decompositions. In this case however the variables are not equiv- 
alent any longer to edge lengths, since as already discussed these do not specify uniquely the 
geometry of polyhedra. Rather, such cellular Regge geometry must use areas and normals as 
fundamental variables. 

Finally, let us make some comments on the coherent states themselves. The discussion so 
far is largely independent of the details of the coherent states on T~Lr- All that is required is 
that they are properly peaked on a point in phase space. The states most commonly used are 
the heat-kernel ones of Thiemann and collaborators. Notice that these are not written in terms 
of the LS coherent intertwiners (41). Nevertheless, it was shown in [41] that they do reproduce 
coherent intertwiners in the large area limit. Alternative coherent states based directly on 
coherent intertwiners appear in [42] . These results show that coherent intertwiners can be used 
as building blocks of coherent spin networks. 

5 On the volume operator 

At the classical level, the volume of a polyhedron is a well-defined quantity. In this section we 
investigate the quantization of this quantity and its relation with the volume operators used in 
loop quantum gravity. 

5.1 The volume of a quantum polyhedron 

Let us consider the phase space Sf of polyhedra with F faces of given area. The volume of the 
polyhedron is a well-defined function on this phase spase, as discussed in Section 3.2. Coherent 
intertwiners provide a natural tool to promote this quantity to an operator in Hp. 

In the following we use the parametrization of the phase space Sf in terms of the cross ratios 
Zk- In particular, the F normals rij are understood as functions of the cross-ratios, rii(Zk). 
Accordingly we call V(ji, Z^) the volume of a polyhedron with faces of area Ai(ji) = 8irjLpji 
and normals rii(Zk), 



For simplicity we assume an ordering of operators such that the area is linear in the spin, 
but the above expression, and the following construction, can be immediately applied to other 
possibilities. 

Let us consider now the Hilbert space of intertwiners Hp associated to the phase space Sf- 
The volume of a quantum polyhedron can be defined in terms of coherent intertwiners Zk) 




(50) 
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and of the classical volume as follows: 

V = j dfj,(Z k ) V(ji, Z k ) \ji, Z k )(ji, Z k \ . (51) 

This integral representation of the operator in terms of its classical version 14 is of the kind 
considered originally by Glauber [45] and Sudarshan [46]. It has a number of interesting 
properties that we now discuss: 

Ql. The operator V is positive semi-definite, i.e. 

(i>\V\iP) = [ dfi(Z k ) V(ji, Z k ) Z k ^)\ 2 > , (52) 



for every \ip) in H. This is a straightforward consequence of the fact that the classical 
volume is a positive function, V(ji, Z k ) > 0. Furthermore, V vanishes for F = 2 and 3. 

Q2. V is a bounded operator in Tip. Its norm 1 1 1 1 = swp^{ip\V\ip) / (ip\ip) is bounded from 
above by the maximum value of the classical volume of a polyhedron with fixed areas, 



dfi(Z k ) V(ji,Z k ) \(ji,Z k \ifj)\ 2 < swp{V(ji,Z k )} = V max (ji) . (53) 
J z h 



Q3. 0-spin consistency. Let us consider the operator V defined on the Hilbert space %f+i 
associated to spins ji, . . . ,jF,jF+i, and the one defined on the Hilbert space Hp associ- 
ated to spins j±, . . . ,jp. When the spin jF+i vanishes, the two operators coincide. This 
is a consequence of the fact that the classical volume of a polyhedron with F + 1 faces 
coincides with the volume of a polyhedron with F faces and the same normals when one 
of the areas is sent to zero. 

These three properties are the quantum version of CI, C2, C3 discussed in Section 3.2. More- 
over, using the fact that for large spins two coherent intertwiners become orthogonal, 

\(j u Z k \j u Z' k )\ 2 ^5(Z k ,Z' k ), (54) 

we have that the expectation value (V) of the volume operator on a coherent state \ji,Z k ) 
reproduces the volume of the classical polyhedron with shape (ji, Z k ), 

(V) = ^f^ff « ViAiWMZk))- (55) 
{Ji, Z k \ji, Z k ) 



14 In the literature [29], the classical function V(ji, Z\) is called the P-symbol of the operator V. On the other 
hand, the expectation value of the operator V on a set of coherent states, i.e. 

Q(ji,Z k ) = {j^ZklV^Zk) , 

is called the Q-symbol. When the P-symbol and the Q-symbol of an operator exist, then the operator is fully 
determined by either of them. The properties of these symbols and of the operator they define have been studied 
by Berezin in [43, 44] 
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This fact allows to estimate the largest eigenvalue of the volume: in the large spin limit, the 
largest eigenvalue is given by V max (Ai), the volume of the largest polyhedron in 5p. 

The spectrum of the operator V can be computed numerically. Let us focus on the case 
F = 4 for concreteness. The matrix elements of V in the conventional recoupling basis are 
given by 

V kk > = {ji,k\V\j h k') = J dfi(Z) V(ji, Z) (j h k\j u Z){ 3i , Z\jt, k'). (56) 

The matrix V k y can be diagonalized numerically to obtain its eigenvalues. 15 We focused for 
simplicity on the simplest case where all the four spins jj are equal to jo- The results using 
Wolfram's Mathematica are shown in Fig. 9 and confirm that the maximum eigenvalue is below 
the volume of the regular tetrahedron. Notice also that the spectrum has a gap. One of the 
interesting questions to investigate in the future is whether this gap survives at higher valence, 
or it decays as for the standard volume operator [47]. 



2.0 - 



1.5 



1.0 
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Figure 9: Some eigenvalues of V. For comparison, the curve is the classical volume of an equilateral 
tetrahedron as a function of the area A = j (units &ir^L 2 p = 1). The empty circles are single eigenvalues, 
the full circles have double degeneracy. The spectrum is gapped and bounded from the above by the 
classical maximal volume, which provides a large spin asymptote. 



It is interesting to notice that the volume operator introduced above commutes with the 
parity operator. This is the operator that sends the normals to their opposite, 

r\j,n) = \j,-n). (57) 

In terms of the stereographic projection, the maps n i— > —n amounts to z i— >■ — 1/z, thus its 
action on coherent intertwiners labeled by the single cross-ratios Z is simply 

V\j i ,Z) = \j i ,Z). (58) 

15 The overlaps, (j,k\ji,Z) = 2^+^+1 [2j+k)\(2j-k)\ -^ / fc( 1 — where is the fc-th Legendre polynomial, and 
the measure, can be found in [9]. 
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Notice that V(Ji, Z k ) = V(Ji, Z k ) thanks to the invariance of the classical volume under parity. 
Moreover the measure d/i(Z k ) is invariant under the transformation Z k — > Z k . As a result, the 
operator (51) commutes with parity, 



vvp = J d f i(z)v(j l ,z)\j i ,z)(j l ,z\ = J 

This explains the degeneracies seen in the spectrum. 



dfi(Z)V(ji,Z)\j u Z)(j u Z\=V. (59) 



Clearly, there are other possibilities for the volume of a quantum polyhedron. All of them 
share the same classical limit, but can have a different spectrum for small eigenvalues. An 

interesting variant is V = y \U\, where U is the oriented- volume square operator, defined as 




dfi(Zk) s(Z k )V 2 {ji,Z k ) \ji,Z k ){ji,Z k \. 



(60) 



Here s(Z k ) is the parity of the polyhedron, i.e. s(Z k ) = ±1 and s(Z k ) = —s(Z k ). 

The operator U anticommutes with the parity, and so does V. Therefore, under the as- 
sumption that the spectrum is non-degenerate, we have that the eigenvalues appear in pairs 
±it. In particular, a zero eigenvalue is present when the Hilbert space T~Lf is odd-dimensional. 
This operator is similar in spirit to the the volume of a quantum tetrahedron introduced by 

Barbieri [3], V B = (87r 7 )iL|,^^| Ji • (J 2 x J 3 )|. In Fig. 10 we show some eigenvalues of V 
and a comparison with Vg. 

For more on semiclassical aspects of the spectrum of the volume, see [48]. 
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Figure 10: Left panel. Some eigenvalues of V. For comparison, the curve is the classical volume of an 
equilateral tetrahedron as a function of the area A = j (units ^tt^L 2 p = 1). All but the zero eigenvalue 
have double degeneracies. Right panel. Same region of the spectrum for Barbicri's operator Vb- Notice 
that here the asymptotic curve is the equilateral volume with areas A = y/jjj + 1). 



5.2 LQG volume operator and the quantum polyhedron 

In LQG, the operator associated to the volume of a region in space is a well studied quantity 
[49, 50, 51]. It is defined on the graph Hilbert space T~Lr as a sum over contributions V n from 
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each node n of the graph within the region R, 

Vr(R) = E ^" ( 61 ) 

nCR 

In order to admit a lifting from %r to the full Hilbert space of LQG, the operator Vr(R) has to 
satisfy a number of consistency conditions that go under the name of "cylindrical consistency" 
[52]. In particular, these conditions are satisfied by the operator V n if (i) it commutes with the 
area of dual surfaces, so that V n reduces to an operator on the intertwiner space %^( n ), and (ii) 
it satisfies a 0-spin consistency condition so that the operators defined on different intertwiner 
spaces coincide when these spaces are identified. 

In the previous section we have introduced an operator V n , given by (51) for the given node, 
that satisfies these conditions. Condition (z) holds because by construction the operator acts 
within %F(n)i anci condition [ii) follows from property Q3 in Section 5. This operator is based 
on the knowledge of the classical system behind the intertwiner space %F(n)- The single node 
operator V n measures the volume of a quantum polyhedron dual to the node, and the operator 
Vr(R) built as in (61) the volume of a region in a twisted geometry. It has a good semiclassical 
limit by construction. 

The standard strategy in LQG is on the other hand rather different. The starting point is 
the classical expression for the volume of a region, 

V(R) = J d 3 x ^j^\ e ijk eabcE a E b E c\ j (62) 

Ef(x) being the Ashtekar-Barbero triad. The key step is to rewrite this quantity in terms of 
fluxes, which are the fundamental operators of the theory. This step introduces a regularization 
procedure which is adapted to a graph T embedded in space. Then, the regularized quantity 
is promoted to an operator in the Hilbert space H-p and the limit of vanishing regulator exists 
and it is well-defined. Two volume operators have been constructed in this way, one by Rovelli- 
Smolin [49], and one by Ashtekar-Lewandowski [50]. Both these operators have the form (61), 
and differ in the regularization procedure and in details on the exact form of V n . For the 
Ashtekar-Lewandowski volume operator, the node contribution is defined on the intertwiner 
space %f as 

V n AL = (STrjf^Lp H V e( ei , ej ,e k ) X ■ {J 3 A J k ) , (63) 




where e(ej,e,-,efe) = ±1,0 is the orientation of the tangents e, to the links at the node. The 
overall coefficient is fixed by a consistency requirement known as 'triad test' [53]. There is a 
large amount of analytical and numerical results on the spectrum of this operator (e.g. [51, 47]), 
particularly because it enters Thiemann's construction of the Hamiltonian constraint [54] and 
thus it is relevant to understand the quantum dynamics of the theory. Moreover its semiclassical 
behaviour has been investigated in detail with the conclusion that only cubulations, that is 
regular graphs with 6-valent nodes, have a good semiclassical limit [55]. In the light of the 
quantum polyhedron introduced in this paper, this result can be understood as follows. 
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On semiclassical states, 16 (Jj) = A4 = Ajrij (see discussion in Section 4 and cf. (37) and 
(42)), and the expectation value of (63) is - at zero order in H [55] 



(Vn L ) = ,h\ E ^ejie^Ai-iljAl^. (64) 

V l<'t<j<fc<F 

As discussed earlier, the variables Ai of the semiclassical state define a polyhedron around the 
node n. The key observation is that (64) is not the volume of that polyhedron. The volume 
of a convex polyhedron with F faces is in general a rather complicated function of the areas 
and normals (see the discussion in Section 3.2). There is however a case where this expression 
simplifies greatly, and in this case it coincides with (64): it happens for parallelepipeds. Paral- 
lelepipeds are a subset of the phase space Sf for F = 6 with areas that are equal in pairs. They 
live within the combinatorial class of cuboids: they are cuboids with three couples of parallel 
faces. 17 The volume of a parallelepiped is 



V=yJ\A 1 -(A 2 /\A 3 )\, (65) 

where (123) are any three faces sharing a vertex. It is straightforward to see that this coincides 
with (64) for the semiclassical state of a cubic analytic node 18 with areas equals in pairs and 
normals parallel pairwise. 

This fact explains why the expectation value of the operator (63) on a semiclassical states 
reproduces the volume of a parallelepiped for F = 6, but not the volume of other polyhedra. 19 



6 On dynamics and spin foams 

Spin foam models for the dynamics of loop quantum gravity are usually built starting from 
a discretization of the spacetime manifold in terms of a simplicial triangulation A. A certain 
control over the dynamics comes from a connection with Regge calculus in the large spin 
limit. Specifically, in this limit the transition amplitudes are related to exponentials of the 
Regge action [9, 31, 56, 57]. This result is generally regarded as a promising step towards 
understanding the low-energy physics of the theory, since discrete general relativity on A is 
reproduced. On the other hand, complete transition amplitudes for LQG require the use of 
more general 2-complexes than those those dual to simplicial manifolds. 20 

16 The semiclassical states used in the analysis of [55] are the heat-kernel coherent states developed by Thiemann 
and collaborators [35]. However, the details on the coherent states do not matter for our argument, all that is 
required is that they are peaked on a given point in the classical phase space Sr. 

17 Notice that parallelepipeds are a set of measure zero among the cuboids. Moreover, cuboids are not the only 
dominant class in phase space Sf with F = 6. 

18 That is, the link are the analytic continuations of each other across the nodes. 

19 It goes without saying that the dependence on areas and normals of the expression (63) can be used to define 
the volume of a tetrahedron, as we saw with Vb earlier. But that would require a different numerical coefficient 
in (63) - an extra y/2/3 - which is hard to motivate in the standard LQG construction. 

20 Although a direct construction of the path integral for arbitrary graph has not been attempted so far, in 
[58] a model valid for arbitrary graph was proposed, based on a natural extension of some algebraic properties 
of the EPR model [59]. 
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Just as Regge calculus is useful to study the semiclassical behaviour on simplicial manifolds, 
a generalization thereof to arbitrary cellular decompositions could be relevant to the full theory, 
and allow us to test whether models such as the one proposed in [58] can be related to (discrete) 
general relativity. In this final Section, we would like to make two remarks on this idea. 

The first remark concerns Regge calculus on arbitrary cellular decompositions. The point 
is that edge lengths are not good variables to capture the (discrete) metric of the manifold. 
This is simply because a generic 4d polyhedron at fixed edge lengths is not rigid. Therefore 
a piecewise-linear metric can not be described by the edge lengths of the polyhedra alone. 
The solution to this problem can be found looking again at Minkowski's theorem, which holds 
in any dimension. The theorem implies that a generic polyhedron in IR™, sometimes called 
an n-polytope, is uniquely characterized by nF — n{n + l)/2 numbers: the volumes of the F 
"faces" (which are now (n — l)-polytopes) and the normals satisfying the n-dimensional closure 
condition. On the other hand, n-simplexes are polytopes with a minimal number of faces, 
F = n + 1. In this case, assigning their n(n + l)/2 edge lengths suffices, thus edge lengths fix 
a unique flat metric on each n-simplex and can be used as fundamental variables in the full 
triangulation. 

Let us fix n = 4. To identify the geometry of each 4-polytope, we need volumes V m and 4d 
unit normals N m of each polyhedron m in its boundary, satisfying the closure condition. For 
these to extend to a piecewise-linear, continuous metric on the whole cellular decomposition, 
we additionally need shape matching conditions, of the sort described in Section 3.4 for three 
dimensions. A tentative Regge-like action can then be written as 

S[V m , N m ] = y~]Af(V m i N m )e f (V m , N m ) + constraints, (66) 
/ 

where / are the 2d faces of the cellular decomposition, and e the deficit angles, defined as usual 
as 2tt minus the sum of dihedral angles of each 4-polytope sharing the face. The constraints are 
the closure and shape matching conditions. In principle, we can interpret (66) as an "effective" 
Regge action in which the internal edge lengths of an initial simplicial triangulation have been 
evaluated on the flat solution. 

The second remark concerns the link between spin foam amplitudes and Regge calculus. A 
lesson from the recent asymptotics studies of the EPR model is that the amplitude is dominated 
by exponentials of the Regge action when the boundary data satisfy certain conditions, which 
guarantee the existence of a unique 4-simplex in the bulk. This suggest that the dominant 
contributions to models on arbitrary graphs could come from requesting the existence of a 
unique 4-polytope, and that the amplitude could be related to a form of the Regge action 
specialized to the 4-polytope, such as the one described above. So the question is whether, 
as for the 4-simplex, the conditions for the existence of the 4-polytope can be mapped into 
conditions on the boundary data, such as 3d closure and non-degeneracy conditions, and shape 
matching. This is a key question that we leave open for future work. We believe that the answer, 
and these considerations in general, will be relevant to tackle the problem of the semiclassical 
limit of spin foams on arbitrary graphs, such as the one proposed in [58]. 
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7 Conclusions 



In this paper we discussed a number of properties of classical polyhedra which are of interest to 
loop quantum gravity. A polyhedron can be uniquely identified by the areas and the normals to 
its faces (Minkowski's theorem [7], Section 2). The identification includes the knowledge of its 
geometry (edge lengths, volume), and its combinatorial class (the adjacency of the faces). This 
information can be explicitly derived from the areas and normals through the reconstruction 
procedure presented in section 3. We observed that the space of polyhedra of given areas is a 
phase space, previously introduced by Kapovich and Millson [6], and used our reconstruction 
algorithm to divide this space into regions corresponding to different classes. 

We then discussed the relevance of polyhedra to the quantum theory. We first recalled 
that the quantization of Kapovich and Millson phase space gives the SU(2)-invariant space of 
intertwiners (section 4), and thus observed that the LS coherent intertwiners can be interpreted 
as semiclassical polyhedra. The polyhedral picture can be extended to a whole graph using 
the twisted-geometry parametrization of the holonomy-flux variables introduced in [1]. The 
knowledge of the classical space behind intertwiners was then used to introduce a new operator, 
which measures the volume of a quantum polyhedron (section 5), and by construction has the 
correct semiclassical limit. We performed some numerical analysis of its spectrum for the 
simplest 4-valent case. We discussed its relation to the volume operators commonly used in 
loop quantum gravity. Finally (section 6), we used the four-dimensional version of Minkowski's 
theorem to make some remarks on Regge calculus on non-simplicial discretizations and its 
possible relevance to spin foam models on graphs of arbitrary valence. 

Our hope is that the notion of a quantum polyhedron can find useful applications in future 
developments of loop quantum gravity, and that the results in this paper are a first step in that 
direction. 
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